Simulating Z2 topological insulators with cold atoms in a one-dimensional optical lattice 



Feng Mei/'^ Shi-Liang Zhu,^ Zhi-Ming Zhang,^'* C. H. Oh,^'"^ and N. Goldman^'^ 

^Laboratory of Photonic Information Technology, LQIT & SIPSE, 
South China Normal University, Guangzhou 510006, China 
^Centre for Quantum Technologies and Department of Physics, 
National University of Singapore, 3 Science Drive 2, Singapore 117543, Singapore 
^Laboratory of Quantum Information Technology and SPTE, South China Normal University, Guangzhou, China 
^Center for Nonlinear Phenomena and Complex Systems - Univer site Libre de Bruxelles , 231, Campus Plaine, B-1050 Brussels, Belgium 

(Dated: January 18, 2012) 

We propose an experimental scheme to simulate and detect the properties of time-reversal invariant topo- 
logical insulators, using cold atoms trapped in one-dimensional bichromatic optical lattices. This system is 
described by a one-dimensional Aubry-Andre model with an additional SU(2) gauge structure, which captures 
the essential properties of a two-dimensional Z2 topological insulator. We demonstrate that topologically pro- 
tected edge states, with opposite spin orientations, can be pumped across the lattice by sweeping a laser phase 
adiabatically. This process constitutes an elegant way to transfer topologically protected quantum states in a 
highly controllable environment. We discuss how density measurements could provide clear signatures of the 
topological phases emanating from our one-dimensional system. 



PACS numbers: 37.10.Jk, 73.43.-f, 67.85.Lm 

The quantum Hall (QH) effect, discovered in 1980, pro- 
vided the first example of a quantum phase that has no sponta- 
neously broken symmetry. Besides, its universal character and 
remarkable robustness have been shown to be related to the 
existence of topological invariants [1-3]. The recent discov- 
ery of the quantum spin Hall (QSH) effect, in materials dis- 
playing strong spin-orbit coupling, has opened the path for a 
new family of topological states: the Z2 topological insulators 
[1,2]. Since then, the search for topological phases of matter 
has become a forefront topic in condensed matter physics [3]. 
In general, topological insulators are insulating in the bulk but 
they feature gapless edge or surface states at their boundary. 
These edge modes are very robust and therefore persist in the 
presence of impurities. The delicate control over these edge 
modes has attracted considerable interest for the realization of 
quantum spintronic and magnetoelectric devices [3]. Further- 
more, in proximity of superconductors, topological insulators 
lead to non-Abelian excitations that could lead to a new archi- 
tecture for topological quantum computation [4] . 

Nowadays, cold atoms trapped in optical lattices are widely 
recognized as powerful experimental tools to mimic a wide 
range of systems originally stemming from condensed matter 
physics [5, 6]. Recently, many experimental efforts have been 
focused on the experimental realization of synthetic mag- 
netic fields and spin-orbit coupling for ultracold atoms [7-1 1], 
which set the stage for the simulation of topological insulators 
and fractional quantum Hall states. In particular, several pro- 
posals have been suggested to realize the QH and QSH states 
using these technologies [12, 13]. 

The QH and QSH phases are realized in 2D systems sub- 
jected, respectively, to strong magnetic or spin-orbit cou- 
plings. Surprisingly, several properties associated to these 
topological states can already be probed through a one- 
dimensional reduction of these systems. This idea has been 
explored theoretically and experimentally in ID quasicrystals 



[14], which reproduce the Hofstadter- Aubry-Andre QH model 
[15, 16]. This elegant discovery opens the possibility to inves- 
tigate QH physics using one-dimensional optical lattices [17]. 

In this paper, we propose an experimental scheme to simu- 
late and detect Z2 topological states with cold atoms trapped 
in a ID optical lattice. We show that these topological phases 
can be described by a generalized ID Harper equation with an 
additional SU(2) gauge structure, which could be simulated 
with a two-component atomic gas trapped in a ID bichromatic 
optical lattice. By adjusting the corresponding laser configu- 
ration, one is able to probe several properties of Z2 topological 
states. In particular, one could transfer the spin-resolved edge 
states from one edge to the other, and measure these states 
through density measurements. One also discusses the pos- 
sibility to define a Z2 topological invariant in this ID frame- 
work, allowing to distinguish between trivial and non-trivial 
topological states. Finally, we describe how density measure- 
ments could provide an efficient tool to measure these invari- 
ants in the present context. 

Let us start by presenting a specific 2D tight-binding model, 
which has been introduced in Ref. [13] to simulate a Z2 topo- 
logical insulator with two-component fermions in an optical 
square lattice. The corresponding second-quantized Hamilto- 
nian reads 

U — fY J JOx^ I JOy{m) I T 

m,n 

+>^stag^( — l)^<^L,^<^m,w, (1) 

where Cm^n is a two-component (spin 1/2) field operator de- 
fined on the lattice site {x = ma, y = na), a = I is the lat- 
tice spacing, and t is the nearest-neighbor hopping amplitude. 
Here, the spin- 1/2 structure derives from the fact that each 
site hosts atoms in two internal states [11]. The second line 
of Eq. (1) describes an on-site staggered potential with am- 
plitude Astag, along the x direction, which has been introduced 
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to drive transitions between different topological phases. The 
Peierls phases 6^ and 6y{m), which accompany the hopping 
along the x and y directions, are engineered within this tight- 
binding model to simulate the analog of spin-orbit couplings, 
and are expressed in terms of the Pauli matrices Oxj^^- The 
space-dependent operator 6y{m) = InmaGz reproduces the 
effect of the intrinsic spin-orbit coupling [18]: it corresponds 
to opposite "magnetic" fluxes ±a for each spin component 
and generates QSH phases. The constant operator = "^^^J^x 
corresponds to a spin-mixing perturbation, and thus simulates 
a Rashba spin-orbit coupling term [18]. 

The Hamiltonian (1) being translationally invariant along 
the y direction, one can explore its properties by imposing pe- 
riodic boundary conditions along this direction. Considering 
this cylindrical geometry, the single-particle wave function is 
expressed as i/A(m,^) = exp(/^^^)^(m), where ky is the quasi- 
momentum along the periodic coordinate, and where the two- 
component wave function = satisfies a gen- 
eralized ID Harper equation [15, 16] 

E^yniky) = t{e'^^^^mMky)^e-'^^^^m-i{ky)) 

^R(m,ky)^^^(ky). (2) 

We have introduced the onsite 2x2 matrix 

R{m^ky) — 2tdisig{cos{27ram -\- ky) ,cos{27ram — ky)) 

+ Astag(-1)"/, (3) 

where / is the identity matrix. The Harper Eq. (2) there- 
fore describes the dimensional reduction of the initial system 
(1), but still captures the essential properties of its topological 
phases. The energy spectrum E = E{ky) obtained by solving 
Eq. (2) displays several energy bands that describe the bulk, 
but also several gapless states which constitute clear signa- 
tures of QSH topological phases [13,18]. We note that Eq. (2) 
generalizes the spinless Harper equation obtained by Aubry- 
Andre, in their study of Anderson localization [16]. Recently, 
the Aubry- Andre model has been simulated in a ID bichro- 
matic optical lattice [19], which motivates us to propose a 
generalization of their setup to investigate the physics stem- 
ming from our spin- 1/2 Eq. (2). 

Let us first focus on the case 7 = Astag = 0, where the 
spin-mixing perturbation Ox and the staggered potential are 
absent. The realization and implication of these terms will 
be discussed at the end of this work. In this case, the spin- 
1/2 model described by Eq. (2) can be realized by trapping 
a two-component atomic gas in a primary ID optical lattice 
Vi{x) = Vi cos^{kix), with wave number ki. In this configu- 
ration, and for a sufficiently deep potential, the atomic system 
is governed by the tight-binding Hamiltonian 

m 

where Cm is a two-component (spin 1/2) field operator defined 
on the lattice site (x = ma) and m = 1 , . . . , L. Then, the on- 
site term characterized by the 2 x 2 matrix R{m) can be real- 
ized by two weak state-dependent lattices, with wave number 



^2, which act independently on the two atomic states [20]. 
For our purpose, the optical potentials associated to these 
state-dependent lattices have the form V|,;(x) = y^cos^(fex± 
0/2), which can be produced by two counterpropagating laser 
beams with linear polarization vectors forming an angle (j) 
[5, 21, 22]. These additional lattices interfere with the pri- 
mary lattice, and supposing that 14 <C Vi , simply lead to the 
onsite perturbation term [19, 23, 24] 

if 2 = A ^ cl^frn^t COS ( 2;rj3 m + ) + c^^c^; COS (27rj3 m - ) , 

m 

(5) 

where A ~ r <C Vi, j3 = fe/^i. Although the phase (j) could 
be affected by an overall shift of Vi {x) with respect to V^^:^{x), 
this parameter can be monitored through various technics (cf. 
the experimental methods in Ref. [24], but also the studies 
of Refs. [25] on the effects of uncontrolled phases in cold- 
atom Aubry- Andre models). Besides, we note that the topo- 
logical properties described in this work remain constant for 
small variations of the parameter (j) (cf. below). In this con- 
figuration, the single-particle equation associated to the total 
Hamiltonian //tot = //i + //2 reads 

^^^(0)=^(^^+i(0)+^^_i(0))+^(m,0)^^(0), (6) 

where 5'(m,0) = Adiag(cos(27rj3m + 0),cos(27rj8m — (j))). 
Therefore, a direct mapping from this ID Aubry- Andre sys- 
tem to the 2D setup of Ref. [13] is obtained by associ- 
ating the commensurability parameter j3 to the "magnetic" 
flux a, the potential strength A to twice the tunneling am- 
plitude 2t, and the phase (j) to the quasimomentum ky. In other 
words, our spin- 1/2 generalization of the Aubry- Andre model, 
which could be simulated using ID state-dependent lattices 
[19], could already reveal the topological properties emanat- 
ing from Eqs. (l)-(2). It is worth emphasizing that, altough 
a direct mapping exists between the energy spectra E{ky) and 
E{(j)), that respectively correspond to the 2D model of Ref. 
[13] and the present Aubry- Andre-type model, the parameter 
(j) is fixed in the latter experimental scheme. Therefore, the 
energy "bands" depicted by the spectrum E{(j)) involves the 
union of different configurations of the system, obtained by 
continuously varying the parameter (j) . 

Let us investigate the spectral properties of Eq. (6): for 
a fixed value of the parameters (j) and A (in the following 
A = 2^ = 2 and (j) G [0,2;r]), and for a rational value of the 
commensurability parameter ji = p/q, the spectrum splits into 
q continua of states (cf. Fig. 1 (a)). These states are delocal- 
ized and describe the bulk of our ID system. In the example 
illustrated in Fig. 1 (a), one has set /3 = 1/3, which leads to 
three "bulk subbands". Between these continua of bulk states, 
and within certain ranges of the parameter 0, one finds two 
degenerate states with opposite spin, whose amplitudes are 
localized at the two edges of the system (cf. the states in Fig. 
1 (a), which are highlighted by a dot (resp. a star) for = 1 
(resp. (j) = 271 — I)). As depicted in Fig. 1 (b), the spin orien- 
tation of these edge states are opposite at the two edges: when 
= 1, one finds a spin-up (resp. down) state atE ^ —t, which 
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Figure 1: (Color online) (a) Energy spectrum as a function of the 
phase for a ID lattice with L = 3S sites. Here (5 = 1/3 and A = 2t. 
(b) The amplitudes |^mtP (resp. are represented in blue 

(resp red) as a function of the site index m, and correspond to the 
three states highlighted in (a) at = l,7r,27r — 1. (c) The spin- 
resolved particle densities p^(m) and p^(m), depicted in blue and 
red respectively, for ^Permi = —t and (j) = 1, 7r,27r — 1. These den- 
sities have been computed for infinitely sharp boundaries. Note the 
inversion of the spin structure at the opposite edges, as (j) is varied. 

is localized at m = L (resp. m = 1). The opposite situation, 
i.e. spin up (resp. down) at m = 1 (resp. m = L), occurs at the 
same energy by setting (j) =2k—1 (cf. Fig. 1 (b)). 

In the analogous 2D system, these two pairs of states con- 
stitute helical edge states, which is a hallmark of the QSH 
effect [18]: when the Fermi energy is set in the first bulk gap 
^Fermi ~ " ^ ^ach edge is populated by opposite spins travel- 
ing in opposite directions. In contrast, in the present context 
of a ID lattice, the parameter (j) is fixed, and therefore, each 
edge is populated by a single spin species. However, by adi- 
abatically varying the phase (j) between [0,27r], following a 
gapless edge state within the lowest "bulk gap", one drives an 
interesting transition between a non-trivial edge-state config- 
uration (e.g. spin up at m = 1 and spin down at m = L) to the 
opposite configuration (e.g. spin up at m = L and spin down 
at m = 1). While the total charge is conserved at each edge, 
this exquisite process leads to spin pumping. We note that 
the edge states remain localized during the whole process, ex- 
cept at singular points (e.g. = tt in the example presented 
in Fig. 1), where these states connect to the bulk. Therefore, 
the edge states survive for small variations of the parameter 
(j), a fact which is in agreement with the topological argument 
discussed below. 

In standard two-dimensional QSH systems, the existence 
of gapless helical edge states inside a bulk gap is guaranteed 
by a topological invariant, the so-called Z2 index (v = 0, 1), 
which can be evaluated from the bulk states [18]. When the 
Fermi energy is fixed in a bulk gap characterized by the non- 



trivial value V = 1, an odd number of helical edge state pairs 
are located at each edge, in which case the system realizes 
the QSH effect. As long as the bulk gap remains open, the 
index v remains constant, which guarantees the robustness of 
the edge states against external perturbations. Now, if one 
maps the example presented in Fig. 1 (a)-(b) to the analogous 
2D system (i.e. (j) ^ ky, where ky takes all the values within 
the range ky G [0, 2;r]), one observes that each bulk gap of the 
spectrum E{ky) hosts two pairs of helical edge states (i.e. one 
pair for each edge). Therefore, the bulk gaps presented in 
Fig. 1 would both correspond to the Z2 index v = 1 in a 2D 
realization of the system [13]. Mapping this example back to 
our ID context (i.e. ky ^ (j), with (j) fixed), this result indicates 
that, as long as the bulk gap remains open, there will always 
be a range G [0i , ^2] between which such edge states will be 
detected. 

Obviously, this topological argument is based on the 
analogy between our ID system and its analogous two- 
dimensional QSH system, whose dimensionality allows to 
properly define the topological invariant V [18]. However, 
one can show that this topological invariant can also be rig- 
orously defined in the ID framework. Indeed, when 7 = 0, 
the Z2 index is simply related to the spin Chem number v = 
SChNmod 2, where SChN = (ChN| - ChN;)/2 and where 
ChN|^^ are the Chern numbers associated to the up and down 
spin respectively [26] . When 7=0, the spin components are 
decoupled and the Chern numbers ChN^^^ can be evaluated 
individually from the standard TKNN expression [27] 

ChNt=;^ £ [dk,d(j>^{\W^;,{k„<t>))), (7) 

— -^Fermi 

where ^ is the Berry curvature associated to the single- 
particle state which is caracterized by the band 
index A situated below the Fermi energy £'Fermi, and where 
kx is the quasimomentum. In this expression, one supposes 
that the parameter evolves continuously along the interval 
[0,2;r]: namely, the definition of the Chern number (7) re- 
quires the union of all the Hamiltonian operators H{^). In 
fact, it was recently shown that such a Chern number could 
be rigorously defined for each [14], and that it remains con- 
stant for all G [0,2;r]. This result, which is in agreement 
with the argument based on the 2D analogy (cf. above), guar- 
antees the existence of edge states for certain ranges of the 
parameter [14]. Finally, we stress that for the general situ- 
ation, where spin-mixing is present 7^0, the topological Z2 
index and the SChN are no longer expressed in terms of the 
individual Chern numbers ChN^^^. However, their values ob- 
tained in the limit 7^0 remain constant for finite 7, as long 
as their associated bulk gap remains open [26]. We note that, 
for the example illustrated in Fig. 1 (a), the Chern numbers 
are given by ChN^^^ = ±1 in the first gap and ChN^^^ = Tl 
in the second, which indeed leads to the non-trivial Z2 index 
V = 1 in both cases. 

One stresses that the presence of spin-polarized edge states, 
for a fixed value of the parameter 0, does not necessarily 
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Figure 2: (Color online) (a) Energy spectrum as a function of the 
phase for a ID lattice with L = 41 sites. Here P = 1/6 and A = 2t. 
(b) The amplitudes l^^tP (rcsp. are represented in blue 

(resp. red) as a function of the site index m, and correspond to the 
four states highlighted in (a), i.e. the states labeled by a dot, a rect- 
angle, a star and a hexagon respectively. 

mean that the Z2 index is non-trivial (v = 1): indeed, one 
should count the number of such states within the whole range 
(j) G [0,27r] to determine whether the number of such helical 
edge state pairs is odd (v = 1) or even (v = 0) at each edge. 
This counting procedure, which could be performed experi- 
mentally by continuously varying the phase 0, allows to rig- 
orously classify the phases of our ID model in terms of the 
Z2 topological index. To illustrate a configuration displaying 
trivial and non-trivial Z2 phases, one has computed the energy 
spectrum and edge state structures for the case j3 = 1 /6 (cf . 
Fig. 2). Here, the first bulk gap is characterized by a non- 
trivial Z2 index v = 1 , as for the previous example discussed 
above for j8 = 1/3. However, in the second bulk gap around 
E ^ —t, one would observe two helical edge state pairs at each 
edge, by varying (j) G [0,2;r] (cf. Fig. 2 (b)): the second gap 
is therefore associated to the trivial phase V = 0. This latter 
result is in agreement with the value SChN = 2, which can be 
computed from Eq. (7) for ^'permi = —t and P = 1/6. 

Now, let us show how the detection of Z2 topological phases 
could already be obtained from density measurements. The 
first strategy would be to directly detect the edge states, which 
is a realistic task for systems featuring infinitely sharp bound- 
aries and low Fermi energy. Indeed, in this configuration, 
the edge states will contribute to the particle density in a de- 
tectable way. We illustrate in Fig. 1 (c) the spin densities, 
defined by 

PuH= I l*Pt,aHI', (8) 

— ^Fermi 

for = 1, TT and 27r — 1 and £^Fermi = —t- The wave-functions 
^t,^,^ (m) are computed from a numerical diagonalization of 
Eq. (6). One clearly observes sharp peaks in the spin densi- 
ties, which correspond to opposite spins at the two edges. As 
(j) is progressively varied from = lto;r = 27r — 1, one first 
observes a transition into the "bulk" regime and finally a spin 
inversion at the edges (cf. Fig. 1 (c)). Let us stress that such a 
result is only valid for abrupt walls: in the presence of an ex- 
ternal confining trap, the edge states and their corresponding 



signatures will be destroyed. In order to overcome this prob- 
lem, one could induce the tunneling t by laser-assisted tunnel- 
ing methods, and create synthetic walls within the confined 
system by abruptly changing the tunneling amplitude in the 
central region [13]. Then, the sharp peaks illustrated in Fig. 
1 (c) and corresponding to the edge states will be observed at 
these synthetic walls. Consequently, by varying the parameter 
0, and performing in situ spin-resolved density measurements 
[28], one could directly detect the spin-pumping process illus- 
trated in Fig. 1 (b)-(c). 

Another method would be to measure the spin Chem num- 
ber SChN, which could also be evaluated from density mea- 
surements. This method is based on the fact that the Chern 
number of individual spin species, ChN^^^, can be computed 
from the density through the Streda formula [29] 

ChN^(^Fermi) = ^, (9) 

Here, the (local) Fermi energy is supposed to lie in a bulk 
gap, which is associated to a plateau in the density profiles, 
and Aj8 = j3 - Ap| = p^iP) - P^iP^)- This equation ex- 
presses the fact that the Chern number associated to a bulk 
gap can be evaluated by comparing the density plateaus ob- 
tained from two configurations of the system (i.e. with j3 and 
j3'). In Fig. 3, one shows the density profiles for j3 = 1/3 
and jS' = 1/4, in a system confined by a harmonic poten- 
tial ytrap(^) = Konf X {m — c)^ , where c = L/2. Consider- 
ing a local-density approximation, one defines the Fermi en- 
ergy locally as ^fermi(^) = ^^fermi " ^trap(w): m this regime, 
where the confining potential is considered to vary smoothly, 
the density profiles depict several plateaus that correspond to 
the bulk gaps located below the chemical potential E^^^^^. In 
Fig. 3, the Fermi energy is set at E^^^^^ = 0, thus the density 
plateaus correspond to the lowest bulk gap of Fig. 1 (a). The 
formula (9) yields ChN|(lst gap) = 1, and therefore v = 1 as 
already evaluated above from Eq. (7). We note that this result 
is independent of the parameter (j), which highlights the fact 
that the topological invariants (i.e. ChN^^^, SChN and v) can 
indeed be defined at each value of the parameter (j) (cf. dis- 
cussion above and Ref. [14]). Consequently, a density mea- 
surement at (j) fixed allows to directly determine the Z2 class 
of our ID system in the presence of an external confining trap. 
Finally, we note that this detection scheme is also suited for 
finite spin-mixing perturbations, i.e. 7 7^ 0, as long as the bulk 
gaps remain open. 

It was shown in Ref. [13], that the combination of a spin- 
mixing perturbation 7 7^ and an additional staggered poten- 
tial Astag 7^ leads to interesting phase transitions between 
trivial and non-trivial Z2 phases. These transitions occur in- 
dividually in the different bulk gaps and can be obtained by 
solving the Harper Eq. (2). Exploring these topological phase 
transitions with our ID model requires to engineer the Peierls 
operator 6^ as well as the staggered potential, which both act 
along the x direction (i.e. the direction of our ID system). The 
realization of the hopping operator Oj demands to control the 
tunneling in a spin-dependent manner, which can be achieved 
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Figure 3: Spin-up density p^(m) /A/3 for a ID lattice with L = 50 
sites. Here j8 = 1/3 (blue line) and = 1/4 (blue dotted line), 
A/3 = /3 — the Fermi energy E^^^^^ = and = 7t/2. The sys- 
tem is confined by a harmonic trap Vconf (m — c)'^, with Vconf = O.Olt. 
Comparing the two plateaus indicates that ChN^ = 1 inside the first 
bulk gap (Eq. (9)): this density measurement attributes a non-trivial 
Z2 index v = 1 to the lowest bulk gap illustrated in Fig. 1 (a). 



with several Raman transitions that act individually on the 
two atomic internal states (cf. Refs. [11, 13]). On the other 
hand, the staggered potential could be easily produced by a 
weak lattice Vstag ~ t, with wave number k^tag = ki/2. Using 
this configuration, one could directly detect topological phase 
transitions by varying the staggered potential strength and per- 
forming spin-resolved density measurements, since the latter 
provide sufficient informations to classify our system in terms 
of the Z2 index v (cf. above). For example, the trivial phase 
V = corresponding to j3 = 1/6, £'Fermi = —t and Astag = 
(cf. Fig. 2) will evolve into a QSH phase with v = 1 for 
Astag>1.5^ (cf. Ref. [13]). 

In summary, we have proposed an experiment scheme to 
simulate Z2 topological phases with cold atoms trapped in a 
ID bichromatic optical lattice. Our scheme is based on the di- 
mensional reduction of a 2D model exhibiting Z2 topological 
insulating phases, and which captures its essential properties. 
Our ID atomic system is described by a generalized Harper 
equation, which can be simulated by a spin- 1/2 generaliza- 
tion of the Aubry-Andre system recently realized with cold 
atoms [19]. The latter can be practically engineered using the 
current technology offered by optical lattices, exploiting the 
interferences of a primary lattice with weak state-dependent 
lattices. Interestingly, our simple scheme is able to trans- 
fer spin-resolved edge states, with opposite spin components, 
from one edge to the other. This manipulation, which can be 
easily performed by varying the secondary lattices configura- 
tion [22], constitutes an elegant manner for transporting topo- 
logically protected quantum states [30]. Besides, we have dis- 
cussed the possibility to drive topological phase transitions, 
by varying an additional staggered potential. Furthermore, we 
have shown that the spin Chern number, which allows to clas- 



sify the topological phases of our ID system, can be eval- 
uated using spin-resolved atomic density measurements. In 
the presence of sharp boundaries, we have shown that spin- 
resolved density profiles would already present clear signa- 
tures of edge states, with opposite spin components at the two 
edges. Therefore, the Z2 phases emanating from our ID sys- 
tem could be probed with the current technologies, such as 
in situ imaging techniques. Let us mention that additional 
signatures could be obtained through cyclotron-Bloch dynam- 
ics [31]. We note that state-dependent lattices generally lead 
to large spontaneous emission rates for fermionic species, a 
drawback which could be avoided by considering an atom- 
chip realization of our model [13]. Finally, we stress that 
our dimensional reduction approach may be generalized to 
explore three-dimensional topological phases [32] using 2D 
optical lattice. 
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